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Abstract 

In this paper, we investigate the effects of random fluctuations of the Earth matter density for long baselines 
on the neutrino oscillation transition probabilities. We especially identify relevant parameters characterizing the 
matter density noise and calculate their effects by averaging over statistical ensembles of a large number of matter 
density profiles. For energies and baselines appropriate to neutrino factories, absolute errors on the relevant 
appearance probabilities are at the level of |APa/ 3 | ~ 10~ (with perhaps lAP^^el/P/ie ^ 1% for neutrinos), 
whereby a modest improvement in understanding of the geophysical data should render such effects unimportant. 
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1. Introduction 

The effects of matter on neutrino oscillations 
|I]-^ in the Earth have been investigated in vari¬ 
ous contexts and with several models |-0- It is 
now well-known that the matter density can sig¬ 
nificantly change the reconstructed neutrino en¬ 
ergy spectrum produced by long baseline neu¬ 
trino experiments, such as by neutrino facto¬ 
ries For most calculations the Pre¬ 

liminary Reference Earth Model (PREM) density 
profile 1^ has been used, which is obtained from 
geophysical seismic wave measurements (see, e.g., 
Refs. [|3^-|^ for information about the structure 
of the Earth’s interior). Furthermore, small er¬ 
rors in the PREM matter density with up to 
5% amplitude have been found and documented 
by many geophysics groups (for a summary, see, 
e.g., Ref. |^). Note, however, that the Earth’s 
matter density distribution is not directly observ¬ 
able from seismological data |3|,^. In this pa¬ 
per, we will discuss how these fluctuations in the 
Earth matter density affect the neutrino oscilla- 
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tion transition probabilities. 

It has been found that for short baselines mat¬ 
ter effects are small Thus, any fluctuations 
could be treated as second order effects and can 
therefore be neglected. For long baselines, how¬ 
ever, the fluctuations may be significant. It was 
noted in Ref. 11 that this effect can be important 
as an additional uncertainty in the determination 
of the CP phase 6cp, especially for certain values 
of 6cp. To estimate this effect, the authors used 
a logarithmic distribution with a certain length 
scale and amplitude with a path integral method 
for the numerical evaluation. On the other hand, 
it was shown in Ref. using a perturbation 
theoretical approach, that fluctuations with small 
amplitudes on length scales much shorter than 
the oscillation length in matter average out and 
give no net effect at all. In this paper, we are in¬ 
terested in the errors on the transition probabili¬ 
ties as functions of the length scale and amplitude 
of the matter density fluctuations. We will, in 
particular, focus on the errors on the appearance 
probability of electron neutrinos and electron an¬ 
tineutrinos at typical neutrino factory energies, 
since matter density noise effects could be rather 
substantial in the determination of the CP phase 
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and matter effects are largest in the appearance 
channel. 

2. The effect of a matter density perturba¬ 
tion 

Before we come to modeling the fluctuations 
in the Earth matter density profile, let us study 
the effect of a single perturbation in the matter 
density. In any quantum mechanical system, de¬ 
scribed by a Schrodinger equation, the impact on 
a free particle’s motion of a potential depends on 
its length scale as well as its amplitude. When 
the length scale of the potential is much shorter 
than the characteristic wave length of the incident 
particle, the particle’s wave function will be un¬ 
able to resolve the exact spatial structure of this 
potential. It is then possible to replace the poten¬ 
tial by a ^-distribution with an amplitude equal to 
the integral of the original potential. Similarly, in 
neutrino oscillations, described by a Schrodinger 
equation, a perturbation on a length scale much 
shorter than the oscillation length in matter could 
be replaced by a ^-distribution. This can be 
seen in Fig. (left plot), where the relative error 
= \Pi.e- P"^f\/PJlf, coming from a 
rectangular matter density perturbation, is plot¬ 
ted as a function of the product of the length scale 
-^/^matter and the density contrast Ap/p. For 
the oscillation parameters we choose 612 = 45°, 
6»23 = 45°, 013 = 5°, Ami2 = 2.5 • 10-3eV^ 
Am^i = 3.65 • 10“^ eV^, and 6 cp = 0, corre¬ 
sponding to the LMA solution with a value for 
013 somewhat below the CHOOZ bound neglect¬ 
ing CP violating effects. The other parameters 
(and terms) are given and described in the figure 
caption. The relative error is shown for several 
fixed values of A <C ^matterj mat¬ 

ter density is varied in such a way that the area 
of the perturbation is unaffected for a fixed value 
on the horizontal axis. For the cut at the verti¬ 
cal line this constant area of the matter density 
perturbation is shown in the right plot, indicat¬ 
ing that the corresponding parameter value is al¬ 
ready a quite pessimistic choice. Apparently, the 
curves in the left plot are approximately equal to 
each other at least below the parameter value in¬ 
dicated by the vertical line and especially for very 


small fixed values of A. This means that the rele¬ 
vant parameter for A is the integral of 

the matter density perturbation, i.e., the product 
of the length scale and the amplitude. Thus, we 
may expect that interference effects, arising from 
non-commuting operators in the Hamiltonian cor¬ 
responding to different matter density layers, be¬ 
come irrelevant for very short length scales com¬ 
pared with the oscillation length in matter. This 
can also be seen in the analytical perturbation 
theoretical approach in App. |^, where the case of 
baselines much shorter than the oscillation length 
leads to first order corrections oc A Ap to the tran¬ 
sition probabilities, i.e., the product of the length 
scale and the amplitude is the relevant parame¬ 
ter. Another result of the numerical analysis is 
that we do not have to take into account isolated 
short scale perturbations, such as the ones com¬ 
ing from the matter density contrast in, for exam¬ 
ple, a mine. Estimating the length scale of such 
a perturbation to be shorter than 10 km and the 
relative density contrast to be of the order of 10%, 
we can read off a relative error much smaller than 
1% from Fig. 0. Since the amplitude of the ap¬ 
pearance probability P^e is basically proportional 
to sin^ 2013 and the error in the determination of 
013 is rather substantial [Q, we will henceforth 
neglect very short isolated perturbations. 

3. A model for the matter density fluctua¬ 
tions 

In this section, we will construct a model for the 
fluctuations in the Earth matter density. Figure ^ 
shows the percentage fluctuations in the Earth 
matter density at a depth of 20 km below the 
Earth’s surface obtained from seismic wave mea¬ 
surements. Though a neutrino beam traverses a 
large range of different depths, one may take this 
figure as an estimate of the characteristic length 
scales and amplitudes involved in the problem. It 
suggests length scales of the order of some thou¬ 
sands of kilometers and amplitudes of the order 
of ±4%, whereas at greater depths one can in¬ 
deed have somewhat larger amplitudes. In ad¬ 
dition, the length scales and amplitudes do not 
seem to vary too much around their average abso¬ 
lute values. This also implies that the transition 
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Figure 1. The plot on the left-hand side shows the relative error in the appearance 

probability of the unperturbed profile with a rectangular matter density perturbation of length 
A and amplitude Ap = p — p, where p = 3 g/cuP is the average matter density. The perturbation 
is assumed to be centered at a baseline of length L = 7400 km. The relative error is plotted as a 
function of the product of the relative length scale and the relative density contrast Ap/p. 

Here T^atter — 17000 km is the oscillation length in matter, determined by the average matter density p 
and the leading neutrino oscillation parameters Am |2 and 0 i 3 for a typical neutrino factory maximum 
energy chosen to he E = 30 GeV. The plot on the right-hand side shows this density perturbation for a 
fixed value of the area spanned by the length scale and the amplitude, z.e., A Ap = 0.01 pL^atten which 
is indicated as a vertical line in the left plot. In both plots, the length scale is fixed to be 5 km (dotted 
curves), 500km (dashed curves), and 2000km (solid curves), respectively. 



Figure 2. Percentage fluctuation in the Earth 
matter density at a depth of 20 km from seis¬ 
mic wave measurements | p4[ (reprinted from 
Ref. |§]). 


regions between negative and positive amplitudes 
are quite short compared to the overall structure. 
So why not simply use these measurements in 
neutrino physics, instead of discussing uncertain¬ 
ties in the Earth matter density? First, these 
measurements contain some averaging as well as 


uncertainties in the equation of state of the Earth 
matter density profile from the seismic wave ve¬ 
locity profile (see, e.g., Refs. I45y46|] ). Second, dif¬ 
ferent groups obtain different results [^, which 
are, however, not qualitatively so much different 
with respect to the parameters we will identify 
below. 

In order to model these fluctuations realisti¬ 
cally and investigate the dependence of the rele¬ 
vant parameters on the neutrino oscillation tran¬ 
sition probabilities, we use a step-function ap¬ 
proach, varying the absolute value of the ampli¬ 
tude Ap > 0 and the length scale A > 0 around 
some average values Apo > 0 and Aq > 0 at ran¬ 
dom. For the random variation we choose a Gaus¬ 
sian distribution with standard deviations ax and 
(TAp, respectively, which is truncated at zero. Fig¬ 
ure 1^ shows some sample profiles for different val¬ 
ues of these standard deviations and for fixed av¬ 
erage values. Comparing this model with Fig. ||, 
we can then estimate a set of realistic parameters 
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Figure 3. Sample profiles for a fixed average length scale Aq and a fixed average amplitude Apo, 
corresponding to the values in Eq. and a baseline length of L = 7400 km. The standard deviations 
are varied relative to the absolute values of the length scale and the amplitude. 


as 


Aq 

2000 km. 

(la) 

ax - 

- 1500 km = 0.75 Ao, 

(lb) 

Apo - 

- 3% p. 

(Ic) 

O'Ap ^ 

. l%p=l/3Apo, 

(Id) 


where p = 3 g/cm^ is the average matter density. 

4. Numerical analysis 

For physical reasons we do not expect surprises 
in certain regions of the parameter space. Since 
we cannot show the results for the whole param¬ 
eter space in Aq, Apo, cr\, and ctap simultane¬ 
ously, we will now systematically investigate the 
dependence on some parameters by keeping the 
other ones fixed. One could do this by show¬ 
ing either the absolute errors in the appearance 


probabilities coming from matter density fluctua¬ 
tions lAP^pl, or the relative errors |APqp|/P^^^. 
For the neutrino channel t'p ^ t'e, the transition 
probabilities are in most regions relatively large. 
In this case, the relative errors are quite mean¬ 
ingful and are usually some percent of the total 
probabilities. Depending on the parameters they 
can sometimes even be larger than 10%. However, 
for the antineutrino channel —>■ Pe the absolute 
probabilities in the denominators of the relative 
errors are rather small. Therefore, it turns out 
that the relative errors are not very sensible in 
this case. Comparing plots for the absolute and 
relative errors and mainly focusing on the quali¬ 
tative behavior, we thus decided to show only the 
absolute error plots. Nevertheless, the fact that 
the absolute values of these errors are rather small 
does not mean that they are small compared to 
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the transition probabilities. For the simulations, 
a large number of matter density profiles was cre¬ 
ated at random and the relative error was aver¬ 
aged over all computations with these profiles. 

Figure ^ shows the absolute errors in the ap¬ 
pearance probabilities for the neutrino channel, 
P^e^ and the antineutrino channel, F/ig, plotted 
as functions of cta for a fixed value of ctap and 
CTAp for a fixed value of cta, respectively. For the 
fixed parameter values we choose, if not otherwise 
noted, the values from Eq. (^. In all of our plots, 
the vertical lines correspond to the parameter val¬ 
ues in Eq. (^, he., the same point in the multi¬ 
dimensional parameter space. Since the curves in 
Fig. ^ are slowly varying close to the vertical lines, 
we will further on take the respective values from 
Eq. (|^) for the standard deviations. In general, 
the absolute errors are growing with larger fluc¬ 
tuations in the length scale and the amplitude. 
This means, in our model, that for more irregu¬ 
lar matter density fluctuations we obtain larger 
errors in the transition probabilities. A possi¬ 
ble explanation could be the zero-truncations of 
the Gaussian distributions, which means that for 
large standard deviations the average values are 
shifted. For antineutrinos the absolute errors are 
in all cases much smaller than for neutrinos, be¬ 
cause for antineutrinos matter effects are, in gen¬ 
eral, much smaller for energies larger than a few 
GeV (no resonance effects). 

Next, let us investigate the dependence of the 
absolute errors on the length scale Aq and the am¬ 
plitude Apo- Here we choose the relative standard 
deviations to be cta = 0.75 Aq and (Tap = 1/3 Apo, 
corresponding to Eq. (|^) for the appropriate Aq 
and Apo. The result of this analysis is shown 
in Fig. H, from which it can be seen that the 
absolute error is essentially proportional to the 
length scale as well as the amplitude of the fluctu¬ 
ations. Indeed, this result is, with respect to the 
average length scale and amplitude, very similar 
to what we obtained from a single perturbation 
in the Earth matter density. In the latter case, 
we observed that the product of the length scale 
and amplitude determines the error in the prob¬ 
abilities. Thus, fixing one of these two parame¬ 
ters gives a linear dependence on the other one. 
This is, of course, only true for a single pertur- 




Figure 4. The absolute errors in the appearance 
probabilities P^e and Ppg, averaged over 20000 
random matter density profiles, with the values 
Aq = 2000 km and Apo = 3% • 3g/cm^ from 
Eq. (|^. The errors are plotted as functions of 
cta/Aq for the fixed value ctap = 1/3 Apo (up¬ 
per plot) and (jAp/Apo for the fixed value a\ = 
0.75 Ao (lower plot), respectively [c/., Eq. (|^]. 
The oscillation parameter values are chosen as de¬ 
scribed in the caption of Fig. ^ 


bation without interference effects, which means 
that corrections to linearity have to be taken into 
account in our more general model. One such 
interesting correction is the bumps in the upper 
plot. Since in our model the number of steps n, in 
which the matter density profile has been divided 
into, depends on the length scale, for a small num¬ 
ber of steps n the (average) transition n —*■ n — 1 
can be seen as a bump in the absolute errors. For 
large n the relative contribution of this effect be- 
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Figure 5. The absolute errors in the appearance 
probabilities and Ppg, averaged over 20000 
random matter density profiles, with the stan¬ 
dard deviation values cta = 0.75 Aq and ctap = 
1/3 Apo, corresponding to Eq. (0)- The errors 
are plotted as functions of the baseline fraction 
Aq jL for three different values of Apo (upper plot) 
and the average matter density fraction Apo / p for 
three different values of Aq (lower plot), respec¬ 
tively. Note that the standard deviations scale 
linearly with the absolute average values. The 
other parameter values are chosen as described 
in the caption of Fig. |^. 


comes negligible. 

Finally, taking the parameter values as given 
in Eq. ( 0 ), the absolute errors in the appearance 
probabilities P^e and Ppg are plotted as functions 
of energy and baseline in Fig. |^. From the upper 




Figure 6. The absolute errors in the appearance 
probabilities P^e and Ppg, averaged over 20000 
random matter density profiles, with the param¬ 
eter values from Eq. The errors are plot¬ 

ted as functions of the energy E for a baseline of 
L = 7400 km (upper plot) and the baseline L for 
an energy of P = 30 GeV (lower plot). The other 
parameter values are chosen as described in the 
caption of Fig. [|. 


plot we can observe that the absolute errors are 
rather small for large energies {E > 10 GeV) and 
so are the relative errors at least for neutrinos. 
Since for a neutrino factory the energy spectrum 
for muon neutrinos is peaked at the maximum en¬ 
ergy 011, the low-energy contributions contain 
less statistical information. This point supports 
our focus on high energies in most of the plots. 

As far as the baseline dependence in the lower 
plot of Fig. H is concerned, the absolute errors 
are vanishing for small baselines, because matter 
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effects become, in general, negligible. Although 
the neutrino and antineutrino curves look quite 
different in this plot, the basic principle of the 
shapes is the same. They are both periodic func¬ 
tions of the baseline, as can easily be seen for the 
antineutrinos, and the difference is only due to 
different period and amplitude. 

5. Summary and conclusions 

In this paper, we hrst investigated the effects of 
a single perturbation in the average matter den¬ 
sity of the Earth (especially the Earth’s mantle). 
We then introduced a model for the description of 
random matter density fluctuations in the Earth’s 
mantle, which is based on the observations ob¬ 
tained from geophysical measurements. Finally, 
we used this model to analyze the dependence of 
the absolute errors in the appearance probabili¬ 
ties Pfj,e and Ppg, especially important for long 
baseline neutrino factory CP measurements, on 
our model parameters. 

We observed that, in particular for neutrinos, 
the absolute error in the appearance probabil¬ 
ity for random matter density fluctuations can 
be quite substantial, corresponding to relative er¬ 
rors of some few percent. The error is essentially 
directly proportional to the product of the ampli¬ 
tude and length scale of a single perturbation or 
the random fluctuations. Furthermore, for short 
baselines L < 1000 km, the errors are vanishing 
for both neutrinos and antineutrinos which means 
that matter density fluctuations can in this case 
be regarded as second order effects compared to 
the small matter effects. 

Finally, we comment on the three possible sce¬ 
narios suggested in the summary of Ref. Q : 


From our analysis, scenario (2) best fits our con¬ 
clusions. In our calculations with randomly gen¬ 
erated matter density prohles, we assumed that 
the fluctuations are completely unknown and we 
obtained relative errors of the order of magnitude 
of some few percent. Comparing the results of the 
measurements of the PREM prohle corrections of 
different groups |^| indicates that there is not 
yet sufficient agreement on the data. Moderate 
reduction of the uncertainty by a more detailed 
analysis of data should help to settle this problem 
for neutrino physics. 


A. A perturbation theoretical approach to 
a single perturbation in the Earth mat¬ 
ter density 

In perturbation theory, we can show that for 
a single perturbation with constant amplitude 
AA = ±2V2GFEAne oc Ap in the average 
Earth matter density A (x p, the first order 
perturbation term is proportional to the “area” 
AS = AAA oc XAp of the perturbation, i.e., the 
length scale A times the amplitude AA. 

The Hamiltonian for the propagation of the 
neutrinos is 

n{r) =nQ+AIC + AAIC = n + AAIC, 

( 2 ) 

where K, = \v^ (z/gl is the projector onto the flavor 
state \ve). Let V(r) = 

be the unperturbed evolution operator, where 
Va = \va){va\ is the projector onto the mass 
eigenstate \va) with the eigenvalue in matter. 
Then the evolution operator of the full evolution 
equation can, to first order in perturbation the¬ 
ory, be written as 


1. The uncertainty of present density models 
poses no signihcant problems. 

2. Moderate reduction of the uncertainty, 
through more detailed analysis of data, is 
required. 

3. Signihcant reduction of this uncertainty, by 
conducting a large scale campaign of geo¬ 
physical observations, is required. 


3 3 

U{r) ~ V(r) - iV{r) AA ^ ^ U Va JC Vb, 

a—1b—1 

(3) 

where 

(4) 

(sa S&) 2 

and ro is the position of the center of the fluc¬ 
tuation. For a huctuation length A, which is 




much shorter than the oscillation length, i.e., 
(^a — ^b)(A/2) 1, the last factor in Eq. d) is 

approximately equal to unity. We then obtain 

(5) 

With this result it is easy to see that the transi¬ 
tion probability Pajsir) at the position r can be 
written as 

PMr) ^ P° 0 {r) + 2ASQ[X*ir)Yir)], 

( 6 ) 

where Papi'^) i® unperturbed transi¬ 

tion probability in constant matter den¬ 
sity, X{r) = {vf)\V{r)\va), and Y(r) = 

{u 0 \ (V(ro -r))'’ We){i^e\V{ro)\va)- This shows 
that the perturbary contribution to the tran¬ 
sition probability is proportional to AS', and 
therefore, largely independent of the form of the 
perturbation. Note that perturbation theory only 
holds for AS small compared to {TC) r ~ Ar, i.e., 
AA/A < r/A. 
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